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ABSTRACT 

New results on the groundstate energy of tight, magnetic knots are presented. Magnetic knots 
are defined as tubular embeddings of the magnetic field in an ideal, perfectly conducting, incom- 
pressible fluid. An orthogonal, curvilinear coordinate system is introduced and the magnetic 
energy is determined by the poloidal and toroidal components of the magnetic field. Standard 
minimization of the magnetic energy is carried out under the usual assumptions of volume- and 
flux-preserving flow, with the additional constraints that the tube cross-section remains circular 
and that the knot length (ropelength) is independent from internal field twist (framing). Under 
these constraints the minimum energy is determined analytically by a new, exact expression, 
function of ropelength and framing. Groundstate energy levels of tight knots are determined 
from ropelength data obtained by the SONO tightening algorithm developed by Pieranski (Pier- 
anski, 1998) and collaborators. Results for torus knots are compared with previous work done 
by Chui & Moffatt (1995), and the groundstate energy spectrum of the first prime knots — up to 
10 crossings — is presented and analyzed in detail. These results demonstrate that ropelength 
and framing determine the spectrum of magnetic knots in tight configuration. 
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1 Introduction 

Work to establish rigorous relationships between energy and topological complexity of physical 
systems is of fundamental importance in both pure and applied mathematics. Progress in this 
direction has been slow, but steady, since Arnold's original contribution of 1974 (Arnold, 1974). 
For magnetic knots, in particular, the problem can be synthesized as follows: if the initial field is 
confined to a single knotted flux tube of signature {V, $} (V magnetic volume and <I> magnetic 
flux), then the minimal magnetic energy M m [ n under a signature-preserving flow is given by 
(Moffatt, 1990) 

M min = m(h)<S> 2 V- l l* , (1) 

where m{h) is a positive dimensionless function of the dimensionless twist parameter h. Of 
particular interest is the value of h for which m(h) is minimal (m m i n ). A fundamental problem 
here is this (Moffatt, 2001): 

Problem 1. Determine m m i n for knots of minimum crossing number 3, 4, 5, ... . 

Minimization of magnetic energy of knot types bears some analogies with another type of 
problem, that originates from work on the shape of ideal knots (see, for instance, the collection 
of papers edited by Stasiak et al., 1998). In this context knots are thought of as made by a 
perfectly flexible and infinitely hard, cylindrical rope closed upon itself; a fundamental question 
here is given by the following (Litherland et al, 1999): 

Problem 2. Can you tie a knot in a one-foot length of one-inch rope? 

This problem admits an obvious generalization to knot types of increasing complexity, so 
that the question can be generalized as to finding the minimal length of a given knot type. As 
we shall see, the two problems tend to coincide at some fundamental level. If the relaxation 
of magnetic field to minimum energy state occurs under a volume- and flux-preserving flow, 
then the process, driven by the action of the Lorentz force, resembles the minimal shortening 
of an infinitely flexible, elastic rope under shrinking deformation. In the incompressible limit, 
a shrinking, volume-preserving flow acts on the tubular knot by increasing the average tube 
cross-section as the knot length diminishes. Thickening of the fattening knot stops when the 
topological barriers given by the knot type prevents further relaxation (see §5 below). For 
magnetic knots this end-state will have minimal energy and, for tight knotted ropes, minimal 
ropelength. Existence of a positive lower bound for magnetic energy (Freedman, 1988), however, 
is not sufficient to guarantee that a global minimum is actually attained, even in ideal conditions. 
With increasing knot complexity, for example, configurational arrangements may indeed prevent 
full relaxation, with local minima of magnetic energy (or ropelength) trapped from further 
minimization. 
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In the present paper we shall consider magnetic relaxation subject to the invariance of 
magnetic signature (volume and flux) with the a priori assumption that the magnetic tubular 
boundary of the flux tube remains circular at all times and that the knot length is independent 
from the internal twist h. There is certainly no physical reason to expect this to happen, other 
than mathematical advantage in the analysis of the energy functional. These assumptions pose 
additional mathematical constraints on the relaxation process, preventing full minimization. 
However, the information on energy thus found provides, we believe, a reasonable approximation 
(from above) to the true bound. By this approach we shall demonstrate (§5.3) that magnetic 
energy minima can be related to the minimal ropelength by an exact, analytical expression 
(given by eq. 39, or the simplified form 54) for the minimized constrained magnetic energy of 
knots. Then, by using minimal length data obtained by the SONO algorithm (briefly reviewed in 
§6.1) developed by Pieranski (1998) and collaborators, we determine the constrained groundstate 
energy levels of the first 250 prime knots up to 10 crossings (§8). In doing so, we also compare 
results extrapolated from the SONO data by using eq. (54) with previous work done by Chui & 
Moffatt (1995; hereafter denoted by CM95), highlighting the limitations of their approach and 
commenting on some marked differences in the results. Some critical issues and open problems 
for future work are discussed in the final section §9. 



2 Magnetic knots as tubular embeddings in ideal fluid 

We consider tubular knots as tubular embeddings of the magnetic field in an ideal, incompress- 
ible, perfectly conducting fluid in S 3 (i.e. IR 3 U {oo}, simply connected). The magnetic field 
B = B(x, t) (x the position vector and t time) is frozen in the fluid and has finite energy, that 
is 

B € {V • B = 0, d t B = V x (u x B), L 2 -norm} . (2) 

A magnetic knot K is given by the embedding of the magnetic field in a regular tubular neigh- 
bourhood T a of radius a > 0, centred on the knot axis C of local radius of curvature p > (see 
Figure 1). The field is actually embedded onto nested tori % (i = 1, . . . , n) in T a , and regularity 
is ensured by taking a < p pointwise along C. The existence of non-self- intersecting nested tori 
in T a is guaranteed by the tubular neighbourhood theorem (Spivak, 1979). C is assumed to 
be a C 2 -smooth, closed loop (sub manifold of 5 3 homeomorphic to S 1 ), simple (i.e. non-self- 
intersecting) and parametrized by arc-length. The total length of C is L = L(C). Evidently 
K has the knot type of C, being either trivial, if C (an unknot) bounds a smoothly embedded 
disk, or essential. For simplicity we take T a = C (g> S given by the product of C with the solid 
circular disk S of area A = it a 2 , taken in the cross-sectional plane to C. The total volume 
V = V(T a ) = na 2 L. We assume the tubular boundary dT a = dT (dropping the suffix) remain a 
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Figure 1: Magnetic knot given by the tubular embedding of the magnetic field centred on the 
knot if 7.7. In ideal conditions knot volume V and magnetic flux <I> are conserved quantities. 



magnetic circular, cylindrical surface at all times, of uniform cross-section all along C; denoting 
by &± the unit normal to &T, we have supp(B) := K, where 

T{K) <-> S 3 , with B • ± = on dT . (3) 

The magnetic flux $ through the cross-section S is given by 

$ = / B • v d 2 x , (4) 

Ja(s) 

where now u is the unit normal to S. In ideal conditions the knot K evolves under the action 
of the group of volume- and flux-preserving diffeomorphisms <p : K — > Magnetic energy 
M{t) and magnetic helicity H(t) are two fundamental physical quantities, defined by: 

M(t) = \( ||B|| 2 d 3 x, (5) 

2 JV(K) 

and 

Hit) = [ A-Bd 3 x , (6) 

JV(K) 

where A is the vector (Coulomb) potential associated with B = V x A. As usual, we take 
V • A = in S 3 . For frozen fields helicity is a conserved quantity (Woltjer, 1958), thus H(t) = 
H = constant. Moreover, it is well-known that helicity admits topological interpretation in terms 
of linking numbers, and for a single magnetic knot we have (Berger & Field, 1984; Moffatt &; 
Ricca, 1992): 
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Theorem 1. Let K be an essential magnetic knot in an ideal fluid. Then 

H = Lk<S> 2 = (Wr + Tw)<S> 2 , (7) 

where Lk denotes the Calugdreanu- White linking number ofC with respect to the framing induced 
by the embedding of the B -lines in T. 

The two geometric quantities Wr and Tw are the writhing number and the twist number: 
Wr is a measure of the average coiling and distortion of C in space and depends only on the 
geometry of C, while Tw measures the winding of the field-lines around C, thus depending on the 
embedding of the B-lines within T. Zero-framing of the field lines denotes zero-linking (Lk = 0) 
of these lines with C, providing a reference measure for helicity calculations (see Appendix). 

3 Curvilinear coordinate system 

It is useful to adopt an orthogonal curvilinear system of coordinates centred on C. Let C be 
parametrized by the equation x = X(s), where s is arc-length, with origin s = at some point 
O € C. Let t(s) = X'(s) be the unit tangent to C, prime denoting arc- length derivative. We 
take C to be inflexion-free, then n(s) and b(s) are respectively the standard unit normal and 
binormal to C, with curvature c = c(s) = l/p(s) and torsion r = r(s) given by the standard 
Frenet-Serret equations, i.e. 

t' = cn , n' = -ct + rb , b' = -rn . (8) 

If c(s) = at some point of C, then C has there an inflexion and n(s) is undefined. Since C 
is assumed to be C 2 -smooth, local inflexional configurations will be resolved by a continuity 
argument from either sides of the inflexion point. Thus, without loss of generality, we shall take 
C to be inflexion- free, that is c(s) > 0, Vs G [0, L], so that the Frenet triad {t, n, b} is everywhere 
well-defined on C. A point Ps5 (see Figure 2) is thus given by 

x = X(s) + r cos #n(s) + r sin $b(s) , (9) 

where (r, is the polar coordinate system in the cross-sectional plane. 

Lemma 2. Let T be a regular tubular neighbourhood of C. Then the system of coordinates 
(r, $r, s), where 

#R = #+ f T(s)ds , (10) 
JO 

provides a zero-twist reference system on C in T. 
Proof. Let us consider a second curve C* , of equation 



X*( S ) = X( S ) + rf 2 (s) , 



(11) 
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Figure 2: Relationship between fixed reference frame and the Frenet frame {t, n, b} for a point 
P in the tube cross-section S. 



parallel transfer of C at a distance r in the normal direction along T2 = ncos??* + bsin?9*; C* 
is a "push-off" of C in the normal direction T2. The pair {C,C*} identifies the ribbon 1Z(C,C*) 
of edges X(s) and X*(s), whose twist provides a measure of the winding of C* around C. For a 
generic infinitesimal displacement of the point P € S (see Figure 2) in space, we have 



dx = t + £ 2 — ^ + ds + T 2 d& + T 3 d£ 3 , 



ds 



ds 



(12) 



where £2, £3 denote coordinates along T 2 and T3 = — nsint?* + bcosi?*, respectively. Note that 
the triple {t,T2,T 3 } is orthogonal. By using (8), we have 



dx 



1 Q ^ 

(l-£ 2 ccostf* +&csintf*)t+ (t + — ) (6T3-6T2 



d S + f 2 d6 + t 3 de 3 . (13) 



Hence, the metric is given by 



dx • dx 



(1 - £ 2 c cos $* + e 3 csinr) 2 + {fi + if) ( T + 



The metric is orthogonal in the following cases: 

dtf 



(i) 



r + — = 
ds 



(h) 6d& - &d6 = . 



(ds) 5 



(14) 



(15) 
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Condition (ii) corresponds to a degenerate system of coordinates, since £2 = £3 + est.; condition 
(i) provides zero twist of the reference system {t, T2, T3} everywhere along C, since (cf. Moffatt 
k Ricca, 1992) 

: ' dtf*(s)\ 

(16) 



Tw = lJc( T{s) + 



ds 



ds . 



Integration of (15)— (i) gives 



'(») = -/ T(s)ds + tf ; 
Jo 



(17) 



without loss of generality we take #0 = 0, so that 1Z(C,C*) spans pointwise in the normal 
direction given by the push-off of C*; the pair {C,C*} is untwisted and provides zero-twist for 
the coordinate system. Hence, the metric is orthogonal and it is given by 



and by taking 



we have 



dx • dx = [1 - c(£ 2 costf* - £3 smr)] 2 (ds) 2 + (d£ 2 ) 2 + (d&) 2 , 



& = r cos -&R , £3 = r sin # R , 



dx-dx= [l-crcos(i? i? + r)] 2 (ds) 2 + (dr) 2 + r 2 (dT? R ) 2 . 
By using (17), the independent coordinate is related to the polar angle 1? by 



(18) 

(19) 
(20) 

(21) 
□ 

The orthogonal system (r, $r,s) was originally found by Mercier (1963). In the orthonormal 
basis {e r ,e^ H ,t} the scale factors h±, I12, /13 are given by 

c^x 

— = h\e r = cos i?(i?r, s)n(s) + sin #(i?r, s)b(s) , 
<9x 



tf K = t?+ I* T(s)ds . 
Jo 



dx 



h2^-d R = — r sin#(i?ij, s)n(s) + r cos#($r, s)b(s) , 



(22) 



— = M = (1 - c(s)r cos 0(0 fl ,s))t(s) . 

^ OS 

By setting k = k(r, tJr, s) = 1 — c(s)r cos #(#r, s), the metric tensor is given by 

/ 1 \ 

(9i,j) = (Mi) ■ {hjej) = r 2 

v fc 2 y 

where ei = e r , e 2 = e# fl , 63 = t, with determinant g = r 2 k 2 , and Jacobian 

J = (/iie r • h 2 e^ R ) x /i 3 t = y/g = rk . 



(23) 
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Figure 3: (a) Toroidal and (b) poloidal flux in T. 



Since we assume that r < c 1 = p, we have J > so that the transformation results well-defined. 
Hence, eq. (9) reduces to 

x = X(s) + re r , with $ = $ R - [ r(s) ds . (24) 

Jo 

4 Magnetic field and flux prescription 

Since dT is a magnetic surface, the magnetic field B is given purely in terms of poloidal and 
toroidal components, that is 

B = B P + B T = B# R e 4R + B s i . (25) 

Moreover, since B is diver genceless, we have 

-d(kB# R ) i d(rB s y 



rk 



+ 



<9??i? ds 
so that there exists a flux function ip(r, i?^, s), such that 



, (26) 



1 dip \ dip 

B *» ~ k~dS ' Bs " • (27) 

For the magnetic field to be single- valued, ip(r, s) must be the sum of terms that are either 
linear or periodic in and s (cf. Bateman 1978, pp. 127-128); hence, periodicity in both 
coordinates may be prescribed. 

We now need to express the magnetic field in terms of toroidal and poloidal flux. Let $t 
be the toroidal flux through a cross-sectional area of T of radius r (see Figure 3- (a)): 

$ T = f B s rM R df . (28) 
Jo Jo 
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From the second of (27) and the periodicity of tp in we have 
hence 

fl d* r (r) , , , . 

fp{r,'&R,s) = — — — — hci(r, s) . (30) 

The poloidal flux $p through a ribbon area of width r (see Figure 3-(b)) is given by 

rL/k pt 

$ P = / / B# R dfkds . (31) 
Jo Jo 

From the first of (27) and the periodicity of ip in s, we have 

(32) 

so that 

iKr,0 R ,s) = ±*^M + C2 (r,t? fl ) . (33) 
Thus, from (30) and (33), we have 

where ip = ip(r,'&R,s) = c\(r,s) + C2(r, $r) is a single-valued function, periodic in ?9r (with 
period 2tt) and in s (with period L). By eqs. (27), we have 

B = I ' I-dT' 2^^rJ + ^ a? "a^J ' (35) 

the total field is given by the sum of an average field, represented by the first term, plus a 
fluctuating field with zero net flux. 



5 Constrained minimization of magnetic energy 
5.1 Standard flux tube 

Let us specify the relation between toroidal and poloidal flux. Let V r = irr 2 L be the partial 
volume of the tubular neighbourhood of radius r; the ratio of the partial to total volume is 
given by V r /V{T) = (r/a) 2 . Now, let f(r/o) be a monotonically increasing function of r/a; for 
example 

/(r/a)= (£) 7 ' (7> ° ); (36) 
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7 = 2 gives the standard ratio of partial to total volume. Denoting by <3? := &t (a) the total 
flux, we have 



where h denotes the magnetic field framing, given by (27r) -1 times the turns of twist required to 
generate the poloidal field from the toroidal field, starting from $p = 0. A direct calculation of 
helicity in terms of fluxes shows that h is indeed the linking number Lk of the embedded field 
(see Appendix). A standard flux tube (cf. CM95) is defined by taking 7 = 2. 

5.2 Topological bounds on energy minima 

Relaxation of magnetic knots under topological constraints has been studied by several authors, 
including Arnold (1974), Moffatt (1990), Freedman (1988), Freedman k He (1991). Various 
bounds on magnetic energy M[t) and relationships between energy minima M m - m and topological 
complexity of knot type were found by these authors. In particular, for zero-framed knots, Ricca 
(2008) has proven that: 

Theorem 3. Let K be a zero-framed, essential magnetic knot, embedded in an ideal fluid. Then, 
we have 



where c m i n is the topological crossing number of K . 

Note: results of Theorem 3 refer to the definition of magnetic energy given by eq. (5), that 
involves the coefficient 1/2. 

By Theorem 1 Lk = yields H = 0, hence for zero-framed knots the lower bound given 
by inequality (38)-(i) simply reduces to state the positiveness of M(t). Equation (38)-(ii), 
however, puts in relation knot topology and magnetic energy minima through crossing number 
information. By comparing eq. (38)-(ii) with (1) we see that c m j n = m(0). Relationship (38)- 
(ii), though, is not exhaustive, since with the exception of the 3- and the 4-crossing knot, all 
other isotopes have several topologically distinct representatives (see section §8 below). Hence, 
inequality (38)- (ii) cannot help to detect knot types of same c m j n on the mere basis of their 
minimum energy level (see the example of Figure 4). For these cases specific other detectors are 
necessary. 

5.3 Constrained relaxation to groundstate energy 

We now look for groundstate energy of magnetic knots by the minimization of magnetic energy 
(5) under specific constraints, following the same approach of CM95. The minimum energy thus 
obtained provides an upper bound to the true minimum attainable in absence of constraints. Let 




(37) 




(38) 
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Figure 4: On the basis of eq. (38)-(ii) two different zero-framed knot types with same volume, 
flux and c m \ n (for example the and K5.2 shown, with c m \ n = 5), may both relax to the same 
minimum energy. 

M* and L* denote the constrained minimum magnetic energy and the corresponding minimal 
knot length. We have. 

Theorem 4. Let K be an essential magnetic knot in an ideal fluid, with signature {V, <&} and 
magnetic field given by (35). We assume that 

(i) {V, <£} is invariant; 

(ii) the circular cross-section S is independent of the arc-length s; 
(Hi) iJj is independent of the arc-length s; 

(iv) the knot length L (or L* ) is independent of the knot framing h. 
Minimization of magnetic energy, constrained by (i)-(iv), yields 



Proof. Let us express the magnetic energy (5) in the orthogonal coordinates (r, i?#,s): 

M(t) = - / / {Bj R + B 2 S ) dr rd$ R kds . 
2 Jo Jo Jo 




(39) 



Since ip is periodic in we define A(r, $r) 



dip/dfiji so to have: 



/ A(r,0 R )d0 R = O 



(41) 
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By assumption (iii) and eq. (41), we have: 



M(t) 



I rL/k r2ir 



JO JO 
1 fL/k i-2-k ra 



1 d$ P \ 2 / 1 d$ T 



,L dr J \2irr dr J 

1 r^/K r'Zir ra 

= - / / (A 2 + 2PA + Q)drrd^Ms , 

2 JO JO JO 



dr rd$R kds 



where 



p 1 (I'l'y ! fd<P r 



i /d* T y 

47r 2 r 2 \ dr / 



(42) 



(43) 



2irr dr ~" L 2 \ dr 

Now we partially minimize M with respect to A, subject to the above constraints; let us introduce 
a Lagrange multiplier A(r) and minimize 

1 



^(A 2 + 2PA + Q)-A(r)A(r,^) 
with respect to A; we obtain 

A = A(r)-P, / (A(r) - P) d# R = 
Jo 

From the second of (45), we have 

1 r 27T 



1 f Z7r 
A(r) = — / PdV R 

27T JO 



Now, by substituting the first of (45) into (42), we get 

^ rL* /k r2ir ra 



1 rLi /K rzir ra . . 

M* = - / / (A 2 (r) -P 2 + Q) drrd^ R Ms , 

2 Jo Jo Jo y ' 



lo Jo Jo 

where (•)* represents constrained minimization; from (43) and (46), we have 

1 r 



o 



L* /d$ T \ 2 2vrr /d$ P x 2 



27rr V dr / L* \ dr 



dr . 



By taking <&p = h$T (cf. eq. 37), the integral (48) reduces to 

M * = ^/d^\ 2 L^+4^W dr _ 
jo V dr / AirrL* 

Finally, by first of (37) and straightforward integration, we have 



M 



+ 



87r(7-l)a 2 2L* I .8(7- 1)F 2L 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



where, in the last term, we have substituted a = (V/irL*) 1 / 2 to express M* in terms of the 
signature {V, $}. 
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□ 

As we see, the constrained minimum energy depends on the length axis L* and the square of 
the twist parameter h. For 7 = 2, eq. (50) reduces further to 

M = 1 $ 2 . (51) 

\2V L* J K J 

At minimum energy, L* is the minimal tube length of the relaxed "tight" knot. 



6 Constrained minimum energy by the SONO algorithm 

As we see from eq. (51), the constrained minimum energy of tubular knots is dictated by the 
tube minimal length. A useful measure of knot complexity is the knot ropelength L/R, where 
from now on (dropping the suffix (•)*) L and R denote the minimal knot length and the radius 
of the maximal, circular cross-section of the tightest knot configuration. For simplicity let us set 
V = 1 and $ = 1. Under signature-preserving flow from V = 1 = itR 2 L, by elementary algebra, 
we have 

^ R ^wm)' 3 ' <52> 



R 



so that we can express the tight knot length L in terms of the ropelength L/R: 

Equation (51) can thus be re-written purely in terms of ropelength and internal twist, that is 

We can now calculate M* by using algorithms that minimize ropelength for any given internal 
twist value. For this we shall rely on one of the most advanced algorithms available, namely the 
SONO algorithm, originally developed by Pieranski (1998) and collaborators. 

6.1 The SONO algorithm 

The Shrink-On-No-Overlaps (SONO) algorithm is a numerical software implemented by Syl- 
wester Przbyl (2001) under the direction of Piotr Pieranski (Pieranski, 1998), and subsequently 
improved by other collaborators. Technical details of this algorithm are widely available in liter- 
ature; the interested reader can consult several papers such as Pieranski et al. (2001), Pieranski 
& Przybyl (2002), Baranska et al. (2004), Baranska et al. (2005). We present here the ideas 
behind the algorithm and we highlight some critical issues. 
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6.1.1 The perfect rope model 

The tightening process is based on the idea that the tubular knot is modelled by a perfectly 
flexible and infinitely hard, cylindrical rope. This means that the rope can be bent with zero 
force (no elastic energy stored) and it cannot be squeezed. This means that the rope cross-section 
remains always perfectly circular. It is also assumed that the rope surface is perfectly frictionless 
and that the knot axis, together with its tangent, is everywhere smooth and continuous. This 
guarantees that the perpendicular cross-section, given by a disk centered on the rope axis, is 
everywhere always well-defined. Moreover, there is a control that enforces that perpendicular 
cross-sections do not overlap to avoid self-intersections. 

6.1.2 Discretization 

The knot axis is standardly discretized piecewise linearly by a finite sequence of segments (beads) 
and vertices, to form a self-avoiding, polygonal knot Kp in IR 3 . The tightening process starts 
from a non-equilateral polygon and proceeds to an equilateral configuration. The main difficulty 
here is to produce a polygonal knot, that best reproduces the knot axis of the perfect rope. In 
finding the tightest (or ideal) configuration the initial ropelength provides an upper bound for 
the minimal value of the ideal configuration. To simulate the hard shell of the perfect rope, 
hard spheres are centred on the Kp vertices. Hence, the collection (union) of all these spheres 
generates a corrugated rope surface. This corrugation is made vanishingly small by increasing 
the number of vertices. 

6.1.3 Original procedure and subsequent improvements 

SONO's basic goal is to minimize the length of Kp, while ensuring that (a) the rope cells 
remain well-defined at all times, and that (b) the rope cells do not overlap, while reducing Kp 
to the equilateral polygon (a rope cell is made by a sphere with two opposite cups appropriately 
removed). To avoid overlapping of contiguous cells a "Control Curvature" (CC) procedure is 
implemented to control (bound) the angle between consecutive beads. This is complemented 
by a second procedure ("Remove Overlaps", RO), that controls the Euclidean distance between 
neighbouring (but not consecutive) vertices. In order to produce an equilateral polygon a third 
procedure ("Equalize Edges", EE) is implemented. Hence, the tightening process is achieved 
by a an iterative application of the EE and RO procedures to minimize the differences in bead 
lengths and take account of emerging overlaps. The tightest configuration is thus given by an 
overlap-free, corrugated rope, centred on an equilateral polygonal knot. A preliminary version 
of SONO (Przybyl, 2001) was followed by an improved version (Baranska et ai, 2004), where a 
finer repositioning of the vertices is implemented through an appropriate displacement towards 
the local center of curvature. 
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6.2 The problem of local/global minimum 

The tightening process forces the rope to change its configuration. The end state is achieved 
when further shortening of the ropelength is no longer possible, because this would create non- 
removable overlaps. When the number of beads is small, the rope surface is highly corrugated, 
and, upon tightening, this could allow different part of the rope to be mutually entrenched in the 
grooves of the corrugation. This would then give false information on the end state, and more 
specifically on local minima of the ropelength (read "energy"). As long as the number of cells 
is high, fine corrugation seems to prevent convergence to unrealistic local minima. However, at 
this stage a further readjustment of the vertices is not allowed and no procedure is implemented 
to explore possible neighbouring minima. Since global minimum configuration is not known a 
priori (and no analytical criterium is available), there is no certainty that the minimal rope- 
length configuration realized by the SONO tightening process corresponds actually to a global 
minimum. 

7 Constrained minimum energy of tight torus knots: compara- 
tive results 

Orthogonal coordinates offer several advantages, one evidently being the economy and trans- 
parency of the calculations involved. Direct comparison of eq. (51) with eqs. (7.1), (7.12) and 
(7.13) of CM95 provides further evidence of this. Before comparing results, however, we should 
point out that eqs. (5.19), (7.4), (7.9), (7.10) and (7.13) of CM95 contains typos and errors that 
need correction: the second integral in the r.h.s. of (7.10), for instance, must be multiplied by 
L to be correct (also on dimensional grounds!). We can then proceed with the comparison: let 
e = a/L; from eq. (51), we have 



By comparing eq. (56) with eq. (7.1) of CM95 (supplied by (7.12) and the corrected (7.13)) we 
can see that the exact expression given by eq. (56) corresponds to the contribution of (7.12) and 
of one term only of the corrected (7.13), with higher-order terms in the series expansion (7.1) 
left out. Here we should point out that the overall contribution of the cut-off terms, being these 
positive and negative, may well sum up to zero. In any case, the result of (7.1) is then applied 




(55) 



which, for 7 = 2, reduces to 




(56) 
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Figure 5: Groundstate energy spectrum of torus knots obtained by Chui & Moffatt (1995). 

to study the groundstate energy of torus knot flux tubes (see §8 of CM95). Here the difficulty 
to analyse tight knot configurations was overcome by using standard torus knot equations (eq. 
(8.1) of CM95) and estimates on contact condition by some approximated form function (eq. 
(8.5) of CM95). The results obtained there are reproduced in Figure 5. 

In considering the relaxation to tight knot configuration, Chui & Moffatt made the implicit 
assumption that throughout the tightening process the knot axis would have not changed its 
initial geometry, the only change in time being an average increase of the tubular knot circular 
cross-section. However, only now we know that this assumption is actually not legitimate. 
Numerical implementation of the global curvature concept to study tight knot configurations 
(Gonzalez & Maddocks, 1999) and extensive numerical simulations made more recently by using 
the SONO algorithm (Pieranski & Przybyl, 2002) show that the knot axis changes shape during 
the tightening process, with progressive deformation due to the change of curvature and torsion 
of C from the initial configuration to the final end-state (while keeping the cross-section circular). 

In the present case we can use the data on minimal ropelength of tight torus knots T Ptq 
(p > q; p, q co-prime integers), obtained by the SONO algorithm simulation (see Przybyl, 
2001; pp. 47-49) for h = 0. The results shown in Figure 6 are given by eq. (54) obtained by 
setting V = 1, $ = 1 (consistently with CM95). The curves are based on the SONO data at 
h = and then extrapolated according to eq. (54). Note that, by doing so, we have assumed 
that minimal ropelength data remain independent from h. The quadratic dependence on the 
twist parameter h is thus evidenced by the family of parabolas. For direct comparison with 
results of CM95 (Figure 5), we report only the constrained energy of torus knots of type T2 tQ 
{q = 3, 5, 7, 9); similar parabolic curves (at higher energy levels) are obtained for torus knots 
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Figure 6: Constrained energy minima of torus knots given by eq. (54). Data obtained by knot 
tightening performed by the SONO algorithm correspond to h = (Przybyl, 2001). 



of higher topological complexity (see §8), such as the T^^ (q > 3), T^ q (g > 4), and so forth. 

The two diagrams show some marked differences and some qualitative similarities. The 
first important difference is in the location of the minima, that, in our case, are all centred at 
h = 0. Consistently with eq. (54), a non-zero framing (given by h ^ 0) yields an increase 
of the energy level, due to the contribution of net twist to the groundstate energy. Moreover, 
take for example the energy curve of the trefoil knot of Figure 5: considering that the 
average writhing number of the tight configuration is given by Wr ~ 3.41 (see, for example, the 
writhing number values calculated by SONO), the value h « 6 (cf. CM95, p. 626) attained at 
the minimum is hard to justify, even by interpreting h in terms of internal twist (since, by eqs. 
(7), and 68 in Appendix, h = if and only if \Tw\ = \Wr\). Moreover, the values at which 
the T3 5 2, topologically equivalent to the Tz^, attains its minimum energy (that with respect 
to T2 ; 3 are different in h and M*) seem to be completely wrong, considering that SONO tests 
on torus knots demonstrated (Pieranski, 1998) that, for given p and q, T q<p and T Ptq relax to 
the same tight configuration (hence, same ropelength and energy). Another difference is in the 
symmetry of the diagrams. The curves shown in Figure 6 for ft / are extrapolated from 
SONO data at h = according to eq. (54). The assumption we made on the independence 
of ropelength from framing prevents us from stating that symmetry is actually preserved when 
this assumption is relaxed. Consistently with the diagrams of CM95, though, the diagrams of 
Figure 6 show that above a certain internal twist threshold h CT = h CT (T pq ), minimum energy 
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Figure 7: Constrained energy minima of prime knots of increasing complexity (c m i n ) given by 
eq. (54). Data obtained by knot tightening performed by the SONO algorithm correspond to 
h = (Przybyl, 2001). Overbar refers to data obtained by averaging values of tight knot lengths 
over the number of representatives within each family. 

levels swap, reflecting the interplay between "internal" and "external" topological complexity 
- the former being associated with the field twist and the latter with the knot embedding in 
S 3 . At zero framing (h = 0), however, internal twist vanishes and energy levels relates to knot 
topology in good agreement with (38)— (ii). 

8 Constrained groundstate energy of prime knots up to 10 cross- 
ings 

We can estimate the constrained energy minima of prime knots up to 10 crossings by using 
data obtained by the SONO tightening process (see Przybyl, 2001; pp. 47-49). In doing so, 
a word of caution is perhaps necessary: as remarked by Przybyl and Pieranski^, one should 
bear in mind that the computed tight knot values tabulated are not rigorous and are subject to 
computational improvements. Assuming that the first decimal digit is correct, we inspect knot 
types of increasing complexity given by their c m i n . For 3 and 4 crossings there is only one knot 
type, but there are 2 for c m j n = 5, 3 for c m \ n = 6, 7 for c m \ a = 7, 21 for c m i n = 8, 49 for c m \ a = 9 

See: http : //f izyka.phys .put .poznan.pl/~pieransk/TablesUpTo9 .html 
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Figure 8: Constrained energy minima of knot representatives (on the right) are compared with 
the corresponding bands of energy (on the left) obtained for h = 0. Results are based on the 
data currently available, obtained by the SONO tightening algorithm at h = (Przybyl, 2001). 
Bands refer to the distribution of energy levels of knot representatives belonging to a given knot 
family, denoted by K Cmin . There is only one knot type for the 3 and the 4 crossing families, 
whereas there are 3 for c m in = 6, 7 for c m in = 7, 21 for c m in = 8, 49 for 

Cmin — 9 and 166 for 

Cmin — 10' 

and 166 for c m in = 10. Since the distribution of the minimal length data produced by SONO is 
not, apparently, a monotonic function of c m in; the corresponding energy levels result sometimes 
swapped between knot types of increasing c m m* To analyse data effectively we proceed in two 
steps. First, from each family of given c m in? we take the average minimal length (L^) out of 
the total number of knot representatives, and use this L to calculate the energy level of that 
family. The resulting curves, shown in Figure 7, are therefore obtained by using eq. (54) and 
extrapolating the diagrams from the SONO minimal average ropelength data obtained at h = 0. 

To analyse the actual distribution of knots in relation to the groundstate energy we present 
a second diagram (see Figure 8), where the median energy levels of Figure 7 are reproduced on 
the right-hand-side for reference. On the left-hand-side we report the actual distribution of 
energy obtained for h = 0, centred on the relative average value (dashed lines). Each band 
shows the spread of energy for each knot family of given c m in (denoted by K Cinin ). For h ^ 
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these bands extend on either side of the vertical axis and are centred on their parabolic median 
(not shown in Figure). 
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Figure 9: Top diagram: constrained groundstate energy of tight knots up to 10 crossings, based 
on data provided by the SONO algorithm; the knots are distributed according to standard 
knot tabulation with increasing crossing number. The table on the right-hand-side lists the 
knots (outliers) with minimum and maximum M*(0) (M* at h = 0) within each c m i n -family. 
Bottom diagram: constrained groundstate energy versus knot type; the knots are now distributed 
according to their minimal ropelength values. 

According to current SONO data, one-to-one correspondence between constrained mini- 
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mum energy and knot complexity holds up to knot K^.is (i- e - it extends to include the first 32 
knot types): the minimal length of Kg^g appears to be shorter than that of K-jx- It is curious 
to notice that, according to standard (Rolfsen) knot tabulation, the knot representatives that 
break this correspondence are those (mostly non-alternating) listed in the last positions of the 
table: these are the last 3 in the 8-crossing family, the last 8 in the 9-crossing family, and, with 
the exception of 3 knots, the last 40 in the 10-crossing family. This is evidenced by the top 
diagram of Figure 9, that shows the energy distribution against knot types listed according to 
the standard knot tabulation. The table on the right-hand-side lists the knot outliers (of each 
c m i n -family) with minimum and maximum M*(h = 0). By re-ordering the knot types of each 
c m i n -family in terms of increasing minimal ropelength we obtain the bottom diagram of Figure 
9, that clearly demonstrates that, besides crossing number, ropelength too is a good detector of 
topological complexity. 

9 Concluding remarks 

The results presented in this paper on the groundstate energy of prime knots up to 10 cross- 
ings are based on a new, exact analytical expression for the constrained minima of magnetic 
energy. This is obtained as an improvement of previous work done by Chui & Moffatt (1995), 
by using appropriate orthogonal coordinates and standard variational methods. A magnetic 
knot is identified with a tubular knot filled by a magnetic field, decomposed in poloidal and 
toroidal components. The tubular knot is given by a tube centred on the knot axis with circular 
cross-section, and the framing of the tubular knot is given by the linking of the field lines with 
the knot axis. By assuming volume and flux invariance and uniform, circular cross-section along 
the tube axis at all times, and independence of knot length from framing, we determine the 
constrained minimum magnetic energy functional (eq. 39), and show that for standard flux- 
tube the minimizer reduces further to a simple expression (eq. 54), function of ropelength and 
framing. Constrained minima are local minima for the magnetic relaxation process and provide 
an approximation from above to the actual minima. By using SONO data for tight knots we 
determine the groundstate energy spectrum of the first prime knots up to 10 crossings, and com- 
pare results for torus knots with previous work done by CM95. Comparison of diagrams show 
some marked differences and some qualitative similarities: contrary to CM95, in the present 
case the constrained groundstate energy reaches its minimum at zero- framing (h = 0), and the 
minimizer results an even function of h. This means that the minimizer for a given framing and 
its negative are simply mirror images. 

Relaxation of basic constraints, numerical improvement of the tightening process and 
improved numerical accuracy can all lead to significant improvements in the energy bounds. 
The assumption of circular cross-section certainly represents the strongest limitation to fur- 
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ther minimization of energy and any relaxation of this constraint is likely to bring appreciable 
improvements on the minima. Evidently, this problem is paired with the difficulty of deter- 
mining the "free volume" of the interstices that survive to the minimization, and since these 
are evidently related to the number of crossings (Buck & Simon, 1999), larger adjustments are 
expected for knots of greater topological complexity. Tight knot data may then be sensitive to 
specific implementation of the tightening process. Apart from the test case of the trefoil knot 
(Gonzalez & Maddocks, 1999), there are no information on how different techniques, based on 
SONO, global curvature radius, or else, compare. Here too, significant improvements of energy 
minima may be achieved when particularly complex topologies are tested. Finally, as work on 
trefoil knot demonstrates (Baranska et al, 2004), computational progress in numerical analysis 
(discretization techniques, error control, etc.) improves accuracy and this, in turn, gets reflected 
in corrections to energy levels. 

From a more general viewpoint, then, a comparison between different techniques would be 
very beneficial. Existence of a global minimum in each isotopy class for minimizers that depend 
solely on ropelength has been established by Cantarella et al. (2002), and an extension of this 
result to include framing is, to the best of our knowledge, yet to be done. Alternative minimizers 
(for some "twisto-elastic" energy) that depend on ropelength and framing are not available, so 
that a direct comparison with results obtained by different methods is not possible. The role 
of knot chirality, for example, that in the present context is simply demanded to ropelength 
minimization, might be detected by minimizers of different type. This calls also for more work 
by numerical simulations of knot tightening, for instance by implementing new procedures not 
only for the fine tuning of local minima, but also for investigating the role of framing, and for 
inspecting symmetry issues in relaxed states. 

In this sense the present work should be considered as an important complement to current 
work on minimum ropelength of tight knots (Litherland et al, 1999; Cantarella et al, 2002) and 
properties of ideal knots (Gonzalez & De La Llave, 2003) and for applications, it has obvious 
implications for all those problems in magnetohydrodynamics that involve estimates on energy- 
complexity relations, especially in astrophysics, solar physics and plasma physics. 
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Appendix 

Magnetic helicity, given by (6), can be written in terms of toroidal and poloidal flux. By using 
the orthogonal coordinates (r, $r,s) we show that eq. (6) reduces to eq. (7) with Lk = h. In 
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order to have <3? = in S S \T and helicity gauge- invariant, we take A = (0, A# R , A s ) such that 

rL/k 
10 

Since B = (0, B$ R , B s ), we must have 
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Prom (35), by equating corresponding components, we have 
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Straightforward integration of the last two equations gives 
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also single-valued and periodic in s and Vr; hence, eqs. (60) become 
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+ C 3 (r) R ,s) , 



(57) 
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where C\(i9r,s) and Cs(i9r,s) are constants of integration: condition (57) is satisfied by setting 
Ci("&r, s) = <frp(a)/L, and continuity of A at r = implies C3 ($r, s) = 0. Thus, we have 
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where <3?p(r) = $p(a) — <I>p(r) is the complementary poloidal flux. Now let 

fL/k r2w rr* 

A Bd°x = / 

*(K) 

the helicity associated with the tubular neighbourhood of radius r* . As in CM95, by using the 
equations above and direct integration, we have 



r fL/k p2n rr* 

H(r*) = A • B d 3 x = / / / (A# B# + A S B S ) dr fdtfp kds , (64) 

Jv*(K) JO JO Jo 



«ir')-l |^p-*$t Uf, (65) 



o 



with total helicity given by 

f a d$ T/ _ f a d$p(f) 



H = H(a)= — -($p(a)-$ P (f))df+ / — f^$ T df, (66) 
Jo dr Jo dr 



hence 

H = 2 I $ T 



r$ T i?£df. (67) 

Jo dr 

By taking d<£p/dr = /i(d$p/dr), magnetic helicity is given by 

H = 2h [ a ^T-^dr = h$ 2 . (68) 
Jo dr 

Hence, eq. (6) reduces to eq. (7), with h = Lk. 
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